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Abstract 

Recently, the stability of certain topological phases of matter under weak perturbations was 
proven. Here, we present a short, alternate proof of the same result. We consider models of 
topological quantum order for which the unperturbed Hamiltonian Hq can be written as a sum 
of local pairwise commuting projectors on a L>-dimensional lattice. We consider a perturbed 
Hamiltonian H = Hq + V involving a generic perturbation V that can be written as a sum of 
short-range bounded-norm interactions. We prove that if the strength of V is below a constant 
threshold value then H has well-defined spectral bands originating from the low-lying eigenvalues 
of Hq. These bands are separated from the rest of the spectrum and from each other by a 
constant gap. The width of the band originating from the smallest eigenvalue of Hq decays 
faster than any power of the lattice size. 

1 Introduction 

Quantum spin Hamiltonians exhibiting topological quantum order (TQO) have a remarkable 
property that their ground state degeneracy cannot be lifted by generic local perturbations [TJ [2] . 
In addition, the spectral gap above the ground state does not close in a presence of such per- 
turbations. This property is in sharp contrast with the behavior of classical spin Hamiltonians 
such as the 2D Ising model for which the ground state degeneracy is unstable in a presence of 
the external magnetic field. Recently the authors of [3] presented a rigorous proof of the gap 
stability for a large class of Hamiltonians describing TQO. The proof was direct, but very long. 
This paper presents a short technical note giving an alternate proof of the same result. We feel 
that both proofs are worth having, since they are complementary. In particular, the present 
proof is much shorter, but less direct. Since we consider the same models and prove the same 
result, many of the definitions in this paper are identical. 
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The interest in proving stability of a spectral gap under perturbations has increased due 
to recent progress in mathematical physics, where a combination of Lieb-Robinson bounds [U 
[SJ E] and the method of quasi-adiabatic continuation [S] with appropriately chosen filter 
functions now provides powerful techniques for studying the properties of gapped local quantum 
Hamiltonians. Perhaps surprisingly, we will use similar techniques to prove the existence of a 
gap. In contrast, to the best of our knowledge, all previous works addressing the gap stability 
problem relied on perturbative methods such as the cluster expansion for the thermal Gibbs 
state [HI [101 [11], Kirkwood-Thomas expansion [12l [13], or the coupled cluster method [IJ1 
Il5j . Our techniques are particularly well suited to handle topologically ordered ground states, 
although they can also be applied to topological trivial states such as ground states of classical 
Hamiltonians (in the classical case stability is possible only for non-degenerate ground states). 
Further, as shown in [31 [7]j the gap stability implies that many of the topological properties of 
the unperturbed system carry over to the perturbed system. 

We consider a system composed of finite-dimensional quantum particles (qudits) occupying 
sites of a D-dimensional lattice A of linear size L. Suppose the unperturbed Hamiltonian Hq 
can be written as a sum of geometrically local pairwise commuting projectors, 

ACA 

such that its ground subspace P is annihilated by every projector, QaP = 0. We impose two 
extra conditions on Hq and the ground subspace P which are responsible for the topological 
order. In [3], it was shown that neither of these conditions by itself is sufficient for the gap 
stability. Let us first state these conditions informally (see Section [2] for formal definitions) : 

TQO-1: The ground subspace P is a quantum code with a macroscopic distance, 

TQO-2: Local ground subspaces are consistent with the global one. 

We consider a perturbation V that can be written as a sum of local interactions 

r>l AeS{r) 

where S(r) is a set of cubes of linear size r and V Ti a is an operator acting on sites of A. We 
assume that the magnitude of interactions decays exponentially for large r, 

max \\V r A\\ < «/e _/ir , 

AeS(r) 

where J, [i > are some constants independent of L. Our main result is the following theorem. 

Theorem 1. There exist constants Jo,ci > depending only on \x and the spatial dimension 
D such that for all J < Jo the spectrum of Hq + V is contained (up to an overall energy shift) 
in the union of intervals Ufc>o^' w he.re k runs over the spectrum of Hq and Ik is the closed 
interval 

h = [fc(l - d J) - 5, fc(l + cjJ) + S], 
for some 5 bounded by J times a quantity decaying faster than any power of L. 
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Using the quasi-adiabatic continuation operators in [23], this super polynomial bound on 5 
can be improved to a bound by an exponential of a polynomial of L, where the polynomial may 
be taken to be arbitrarily close to linear. 

The stability proof presented in [3] involved two main steps: (i) proving stability for a special 
class of block- diagonal perturbations V composed of local interactions V r a preserving the ground 
subspace of Hq, and (ii) reducing generic local perturbations to block-diagonal perturbations. 
The most technical part of the proof was step (ii) which required complicated convergence 
analysis for Hamiltonian flow equations. 

In the present paper we show how to simplify step (ii) significantly using an exact quasi- 
adiabatic continuation technique [TJ, [8] . In Section [3] we define two types of block-diagonal 
perturbations: the ones that preserve the ground subspace of Hq locally ([V^^P] = 0) and 
globally ([V, P] = 0). We prove stability under locally block-diagonal perturbations in Section [5] 
which mostly follows [3] and uses the technique of relatively bounded operators. The reduction 
from globally block-diagonal to locally block-diagonal perturbations is in two steps; in Lemma[7J 
we exploit the idea of [T7] to write a Hamiltonian of a gapped system as a sum of terms 
such that the ground states are (approximate) eigenvectors of each term separately, and then 
we show that such terms can be written as locally block-diagonal perturbations in Section 
Finally, we use an exact quasi-adiabatic continuation technique presented in Section [6] to reduce 
generic local perturbations to globally block-diagonal perturbations in Section [7J 

In contrast to [3] we use certain self-consistent assumptions on the spectral gap of the 
perturbed Hamiltonian H s = Hq + sV. In Section [7] we prove that if the minimal gap A m j n 
along the path < s < 1 is bounded by some constant (say 1/2) then, in fact, A m j n is bounded 
by a larger constant (say 3/4). We use continuity arguments to translate this result into an 
unconditional constant bound on A m j n , see Section [7J The continuity arguments we use rely on 
the fact that the spectrum of Hq + sV is a continuous function of s for any finite sized system; 
however, our theorem gives bounds that are independent of system size. 



2 Hamiltonians describing TQO 

To simplify notation we shall restrict ourselves to the spatial dimension D = 2. A generalization 
to an arbitrary D is straightforward. Let A = x be a two-dimensional square lattice of 
linear size L with periodic boundary conditions. We assume that every site it € A is occupied 
by a finite-dimensional quantum particle (qudit) such that the Hilbert space describing A is a 
tensor product 

H = (g)H u , dimH u = 0(l). (1) 

Let S(r) be a set of all square blocks A C A of size r x r, where r is a positive integer. Note that 
S(r) contains L 2 translations of some elementary square of size r x r for all r < L, S(L) = A, 
and Sir) = for r > L. We shall assume that the unperturbed Hamiltonian Hq involves only 
2x2 interactions (otherwise consider a coarse-grained lattice): 

A€S(2) 
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Here Qa is an interaction that has support on a square A. We also assume that the interactions 
Qa are pairwise commuting projectors, 

Qa = Qa, QaQb = QbQa for all A,B G 5(2). (3) 

Well known examples of Hamiltonians composed of commuting projectors are Levin- Wen mod- 
els [18] and quantum double models [2]. 

We assume that Hq has zero ground state energy, i.e., the ground states of Hq are zero- 
eigenvectors of every projector Qa- Let P and Q be the projectors onto the ground subspace 
and the excited subspace of Hq, that is, 

P= XI V-Qa) and Q = I-P. (4) 

A£S(2) 

For any square B £ 5(r), r > 2 define local versions of P and Q as 

Pb= II ( J -&0 and Q B = I-P B . (5) 

AeS(2) 

Note that Pb and have support on B. 

We shall need two extra properties of Hq related to TQO defined in [3j. We shall assume 
that there exists an integer L* > L a for some constant a > such that one has the following 
properties: 

1. TQO-1: Let A € S(r) be any square of size r < L* . Let Oa be any operator acting on 
A. Then 

PO A P = cP 

for some complex number c. 

2. TQO-2: Let A € S(r) be any square of size r < L* and let B G S(r + 2) be the square 
that contains A and all nearest neighbors of A. Let Oa be any operator acting on A such 
that O a P = 0. Then 

O a Pb = 0. 

The property TQO-1 is often taken as a definition of TQO, see Ref. [H EJ [HI |20] . One can 
also show that L* > d/A where d is the distance of P considered as a codespace of a quantum 
error correcting code, i.e., the smallest integer m such that erasure of any subset of m particles 
can be corrected, see Ref. [21] for details. 

We will need later a corollary of TQO-1 and TQO-2. We use h{A) to denote the square 
containing square A as well as all first, second,..., Z-th neighbors of square A. That is, it is a 
"ball" of distance / around A. 

Corollary 1. Let Oa be any operator supported on a square A. Let C = 02{A) and suppose 
that C has size bounded by L* . Then 

P c O A Pc = cP c , (6) 
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for some constant c. 

Thus, all states tp with Pcip = ip have the same reduced density matrix on A, and in 
particular the reduced density matrix of i/j on square A is equal to the reduced density matrix of 
any ground state on A. 

Finally, we have the useful equality: 

\\O a P\\ = \\0 A P C \\. (7) 

Proof. Let B = b\{A). Note that PbOaPb commutes with the Hamiltonian, and hence com- 
mutes with P. So, P (P b O a Pb) P = (P B O a Pb) P- So by condition TQO-1, {P B A P B ) P = cP 
for some c. Define O a = A — c. Then, (PbO' a Pb) P = 0. So, by condition TQO-2, 
(P B 0' A P B ) Pc = 0. Note that P B 0' A P B commutes with P c . So, P c (P b O' a P b ) P c = 0. Hence, 
PcO' A P c = 0. So, P c O A P c = cP c . 

Now consider the second claim. Let V be any state such that Pcip = V 7 - Let p A {^) be the 
reduced density matrix of ip on A. Then, tr(p A A ) = tr(Pc\ip){ip\PcO A ) = c(0 A )tr(\ip){ip\) = 
c(0 A ), where the constant c(O a ) depends only on O a and not on -ip. Thus, the reduced density 
matrix is the same for all vectors such that Pcip = V'- 

Finally, let ip, * satisfy ||Oa-Pc|| = \0 A ip\ and ||Oa-P|| = lO^ol- Then, we have that 
P c ip = t/j and P^q = ^o, which implies tx(\0 A \ 2 \ijj) = tx(\0 A \ 2 \^Q)(^>Q\), since their reduced 
density matrices agree on A. Hence, ||Oa-P|[ = ||Oa-Pc||> as claimed. 

□ 

3 Local Hamiltonians and the Lieb- Robinson bounds 

Throughout this paper we restrict ourselves to Hamiltonians with a sufficiently fast spatial 
decay of interactions. Any such Hamiltonian V will be specified using a decomposition into 
local interactions, 

V = H £ Vr,A) (8) 
r>l AeS{r) 

where v} A = V r ,A is an operator acting non-trivially only on a square A. We shall often identify 
a Hamiltonian and the corresponding decomposition unless it may lead to confusions. Let us 
define several important classes of Hamiltonians. 

Definition 1. A Hamiltonian V has support near a site u G A if all squares in the decomposition 
Eq. contain u. 

Definition 2. A Hamiltonian V is globally block- diagonal iff it preserves the ground subspace 
P, that is, [V, P] = 0. A Hamiltonian V is locally block- diagonal iff all terms V rA in the 
decomposition Eq. (EJj preserve the ground subspace, that is, [V r)A ,P] = for all r,A. 

Definition 3. A Hamiltonian V has strength J if there exists a function f : Z + — >• [0, 1] 
decaying faster than any power such that 

\\V r)A \\ < Jf(r) for allr>\ for all A G S(r). 
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Our main technical tool will be the following corollary of the Lieb-Robinson bound for 
systems with interactions decaying slower than exponential |24l 120] . 

Lemma 1. Let H be a (time- dependent) Hamiltonian with strength 0{\). Let V be a Hamil- 
tonian with strength J. Let U(t) be the unitary evolution for time t, \t\ < 1, generated by H . 
Define 

V = U(t)VU(tf. 

Then V has strength 0{J). If V has support near a site u then V also has support near u. 

We will also need an analogue of Lemma[T]for an infinite evolution time. It will be applicable 
only to evolution under Hamiltonians H having a finite Lieb-Robinson velocity [31 El El 123] > for 
example, Hamiltonians with exponentially decaying interactions. 

Definition 4. A Hamiltonian V has strength J and decay rate /i iff 

1 1 ^r, A || — J ex P ( — H r ) f or all r > 1 for all A € S(r). 

Lemma 2. Let H be a (time- dependent) Hamiltonian with strength 0(1) and decay rate [i > 0. 
Let V be a Hamiltonian with strength J. Let U(t) be the unitary evolution for time t generated 
by H . Finally, let g(t) be any function decaying faster than any power for large \t\. Define 

/oo 
dtg(t)U{t)VU(t) ] . 
-oo 

Then V has strength 0(J). If V has support near a site u then V also has support near u. 

The superpolynomially decaying function of r bounding the norm of V t ^a depends on the 
superpolynomially decaying function of r bounding the norm of V Tj a as we h as on \i and the 
function g(t). (Remark: In fact, one can show that there exist g(t) that have the properties 
needed later and that decay as an exponential of a power of t, making it possible also to consider 
Hamiltonians H that do not have a decay constant /x > but that instead have sufficiently fast 
stretched exponential decay. We omit this for simplicity.) 



4 Reduction from global to local block-diagonality 

In this section we prove that a certain class of globally block-diagonal perturbations can be 
reduced to locally block-diagonal perturbations with a small error by simply rewriting the 
Hamiltonian in a different form. 

Lemma 3. Consider a Hamiltonian 

H = H + Y J X U , (9) 

where Hq obeys TQO-1,2 and X u is a perturbation with strength J that has support near u. 
Suppose that 

[X u ,P}=0 for all u. 
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Then we can rewrite 

H = H + V + A, (10) 

where V is a locally block- diagonal perturbation with strength 0{J), and |j A|| decays faster than 
any power of L* . 

Proof. Consider some fixed site u and let X = X u . By assumption, X has a decomposition 
X = ^2i r >iX{r), where X{r) = ^AgS(r) Abu ^r,A and the norm ||X(r)|| decays faster than 
any power of r. By adding constants to the different terms X rt A and using TQO-1 we can 
achieve PX rA P = for all r < L* . Define A = PX = PXP = J2 r >L* PX{r)P. Note that 
|| A || decays faster than any power of L* . Define also X' = X — A such that 

X'P = 0. 

We can assume that X' has strength O(J) and its support is near u if we treat A as a single 
interaction on a square of size L. To simplify notation we set X = X' in the rest of the proof. 
Choosing any / < L* and applying Eq. ([7]) to O = Yll=i ^( r ) we arrive at 

l l 

ii x(r)p h+2{u) ii = ii E X ^ r ) p w ^ w xp w + ii E x ^ r ) p \\ ^ J f® ( n ) 

r=l r=l r>l 

for some function / decaying faster than any power of I. 

Consider a given site u. Let B r be the union of all squares of size r that contain u, such 
that X(r) has support on B r . We have B\ C i?2 C . . . C Bm = A for some integer M. Define 
an orthogonal unity decomposition / = Ylm=i by 

Ei = Qbi ^ 

Em = QB m P Bm ^ for 2 <m<M, 
E M +i = Pb m =P- (12) 

Taking into account that Em+iX = XEm+i = we arrive at 

X= Yl E p X(q)E r = YY(j)+^2z(q), (13) 

l<P,Q,r<M j>l q>l 



where 




Y(j)= E E P E X (l)) E r- (14) 



and 



l<p,r<M 
p+r=j 



Z(q)= E E P X ^) E r- (15) 

l<p,r<g+l 



All operators Y(j) and Z(j) are hermitian, annihilate P, and act only on sites in the square 
Bj + \. The norm of Z(q) decays faster than any power of q because of the decay of the norm of 
X(q). 
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We claim that the norm of Y(j) decays faster than any power of j. Indeed, Y(j) is a sum 
over j — 1 different terms corresponding to different choices of p, r. We show that the norm of 
each such term p, r decays fast enough. Assume without loss of generality that p > r. Then, 
ll-^P {^2^=1^^ 2 X(q)\E r \\ < \\E p (^2™^ p ' r ^ 2 X(q) s j ||, which decays faster than any power of 
p by Eq. (fTT|) . Since p > j/2, it decays faster than any power of j. 

Thus, we have decomposed X as a sum of terms which annihilate P, supported on squares 
of increasing radius about site u, with norm decreasing faster than any power of the size of the 
square. This completes the proof. □ 



5 Stability under locally block- diagonal perturbation 

In this section, we define the concept of relatively bounded perturbations (see Chapter IV of |22j) 
and show that the spectral gaps separating low-lying eigenvalues of Hq are stable against such 
perturbations. We then demonstrate that locally block-diagonal perturbations satisfy the rela- 
tive boundness condition which proves the gap stability for locally block-diagonal perturbations. 

Note that our condition of relative boundness is different from the one used in [10J . In 
particular, our condition provides an elementary proof of the gap stability without the need to 
use cluster expansions as was done in |10| . 

Let Hq and W be any Hamiltonians acting on some Hilbert space %. We shall say that W 
is relatively bounded by Hq iff there exist < b < 1 such that 

||W^|| < b \\H i/f\\ for all \i>) G H. (16) 

The following lemma asserts that a relatively bounded perturbation can change eigenvalues 
of Hq at most by a factor 1 ± b. 

Lemma 4. Suppose W is relatively bounded by Hq. Then the spectrum of Hq-\-W is contained 
in the union of intervals [\q(1 — b), Xq(1 + b)] where \q runs over the spectrum of Hq. 

Proof. Indeed, suppose (Hq + W) = X\tp), that is, 

(H Q -\I)\^) = -W\i>). (17) 
The relative boundness then implies ||(-ffo — -^-OV'II — ^ll-^o^ll) that is, 

mH -xi) 2 ^)<b 2 (mo^)- (is) 

Let Hq = Y2x n ^oP\ be the spectral decomposition of Hq. Here the sum runs over the spectrum 
of Hq and P\ is a projector onto the eigenspace with an eigenvalue Ao- Define a probability 
distribution p(Xq) = {ip\P\ \ip} . Substituting it into Eq. (fT8|) one gets 

]T(Ao-A) 2 p(Ao) <]T6 2 A^(A ). (19) 

Ao Ao 
Therefore there exists at least one eigenvalue Ao such that 

(Ao - A) 2 < 6 2 A 2 . (20) 

This is equivalent to Aq(1 — b) < A < Aq(1 + b). □ 
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In the rest of this section we proof stability under locally block-diagonal perturbations. Let 
Hq be a Hamiltonian satisfying TQO-1,2. 

Lemma 5. Let V be a locally block- diagonal perturbation with strength J. Then the spectrum 
of Hq + V is contained in the union of intervals | s J fc>0 Ik, where k runs over the spectrum of Hq 
and 

I k = {A € R : fc(l - b) - 5 < A < k{l + b) + 5} 
for some b = O(J) and for some 5 decaying faster than any power of L* . 

Proof. For any integer r > 1 define 

W(r)= y rA- 

AeS(r) 

By assumptions of the lemma we have 

\\V r ,A\\ < Jf(r) (21) 

for some function / decaying faster than any power and [V^4,P] = 0. Performing an overall 
energy shift and using TQO-1 we can assume that 

V rA P = PV r ,A = for all 1 < r < L* and for all A G S(r) (22) 

Proposition 1. Suppose 1 < r < L* . Then W(r) is relatively bounded by Hq with a constant 
b(r) = 0{Jr 2 f{r)). 

This proposition immediately implies that ]Ci<r<L* ^( r ) is relatively bounded by Hq with a 
constant b = Ylr=i K r ) = 0{J). Lemmad]then implies that the spectrum of -f^o+Z]i<r<L* ^( r ) 
is contained in the union of intervals [k(l — b), k(l + b)]. Treating the residual terms 2r>L* ^( r ) 
using the standard perturbation theory then leads to the desired result. 

Proof of the proposition. Let A = B\ U Bi U . . . U E>m be a partition of the lattice into 
contiguous squares of size r such that any 2x2 square is contained in exactly one square B a . 
If a 2 x 2 square is on the boundary of two, or more, of the B a , we include it in the bottom 
left one in the partition. Moreover, if r does not divide L, then squares at the boundary may 
be truncated to rectangles, to fill in the partition. We refer to B a as boxes to distinguish them 
from the squares involved in the decomposition of V. For any binary string y € {0, 1} M define 
a projector 

M 

Ry=\{[yaQB a + {l-ya)PB a \. 
a=l 

Clearly the family of projectors Ry defines an orthogonal decomposition of the Hilbert space, 
that is, ^2y Ry = I. Given a string 3^ € {0, 1} M , we shall say that a box B a is occupied iff 

y a = i- 



9 



We claim that any operator V V) a acting on a square A 6 S(r) and satisfying Eq. (f22j) has 

only a few off-diagonal blocks with respect to this decomposition. Specifically, TQO-2 implies 
that 

RyV r , A Rz + (23) 

only if A has distance 0(1) from some occupied box in y and A has distance 0(1) from some 
occupied box in Z, and the configurations y, Z differ only at those boxes that overlap with A. 
Clearly, for any fixed y such that y has k occupied boxes the number of pairs (A € S(r),Z) 
that could satisfy Eq. (|23|) is at most 0(kr 2 ). Denote for simplicity 



W = W(r), w = Jf(r). 



For any state 1^) we get 



\W 2 \ip) = (ip\RyWR z WR v \ip) 
y,z,vc[M] 

Y \\RyWRz\\ • \\R Z WR V \\ ■ \\R y ^\\ ■ \\R v ^j\\ 



< 



< 



y,z,v 

WRyWRzW ■ \\RzWR v \\ ■ \{Q>\Ry\il>) + (ip\RvW) 



y,z,v 



= \\RyWRz\\ ■ \\R 2 WR V \\ ■ (i/}\Ry\ip) 

y,z,v 

< E E 0(k 2 r i w 2 )^\Ry^)=0(w 2 r^\G^}, (24) 
fc>o y-.\y\=k 



where 



G = E E k2R y- ( 25 ) 

fc>o y-.\y\=k 

The inequality Eq. (|24p follows from the fact that y and Z differ at at most O(l) boxes and an 
obvious bound k(k + 0(1)) = 0(k 2 ). Finally, note that G < Hq since any configuration with k 
occupied boxes must have at least k defects (squares A € 5(2) for which Qa has eigenvalue 1) 
and since creating a defect costs at least a unit of energy. We arrive at 

(ij\W 2 \if>) < b 2 {^\H 2 \il;), b = 0(wr 2 ). (26) 

This completes the proof. 

□ 
□ 

6 Exact quasi-adiabatic continuation 

We define a continuous family of Hamiltonians, 

H S = H + sV, (27) 
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so that as s varies from to 1, H s continuously interpolates between Hq and the perturbed 
Hamiltonian. 

We define a quasi-adiabatic continuation operator, T> s by 

W s = J dtF(t)exp(iH s t)(d s H s ^jexp(-iH s t), (28) 

where the function F(t) is defined to have the following properties. First, the Fourier transform 
of F(t), which we denote F(lj), obeys 

\u\ > 1/2 -> F(oj) = (29) 

Second, F(uj) is infinitely differentiable, so that F(t) decays faster than any power of \t\. Third, 
F(t) = —F(—t), so that T> s is Hermitian. 
We define a unitary operator U s by 



U s = S' exp(i I ds'V s 
Jo 



(30) 



where the notation S' denotes that the above equation (|30p is an s'-ordered exponential. 

The motivation for defining the above unitary operator is contained in the following lemma. 

Lemma 6. Let H s be a differentiable family of Hamiltonians. Let |^ l (s)) denote eigenstates of 
H s with energies E{. Let E m i n (s) < E max (s) be continuous functions of s. Define a projector 
P(s) onto an eigenspace of H s by 

{s),Em ax (s)] 

P(s)= \^(s))(^(s)\. (31) 

i 

Assume that the space that P(s) projects onto is separated from the rest of the spectrum by a 
gap of at least 1/2 for all s with < s < 1. That is, all eigenvalues of H s are either in the 
interval [E m i n (s), E max (s)], or are separated by at least 1/2 from this interval. Then, for all s 
with < s < 1, we have 

P(s) = U s P(0)Ul (32) 

Proof. This is lemma 7.1 in [3]. It is also well-known in the theory of quasi-adiabatic continua- 
tion. □ 



7 Stability proof 

In this section we prove Theorem [TJ For any s € [0, 1] define 

H S = H + sV. 

Let g be the ground state degeneracy of Hq. Choose E m i n (s) as the smallest eigenvalue of H s and 
let E max (s) be the g-th smallest eigenvalue of H s (taking into account multiplicities). Finally, 
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let A(s) be the spectral gap separating eigenvalues of H s in the interval [E m i n (s), E max (s)] from 
the rest of the spectrum. We shall choose the constants Jq,c\ sufficiently small such that the 
interval Jo is separated from 1^, k > by a gap at least 3/4. Then the theorem implies that the 
interval [E m i n (s), E max (s)\ is contained in Iq for all s G [0, 1] and hence 

A(s) > - for all < s < 1. (33) 

Suppose we have already proved the theorem for the special case when 

A(s) > - for all s € [0,1]. (34) 

In the remaining case there must exist s* E [0,1) such that A(s) > 1/2 for s 6 [0, s*] and 
A(s*) = 1/2 (use the fact that A(0) > 1 and continuity of A(s)). Applying the theorem 
to a perturbation s*V which satisfies Eq. (f3"4"|) we conclude that A(s*) > 3/4 obtaining a 
contradiction. Thus it suffices to prove the theorem for the case Eq. ()34p . 
Define 

H' s = U}(H + sV)U s = U\H S U S , (35) 

where 

f/ s = S' exp(i / ds'X» s /), (36) 
Jo 

is the exact adiabatic continuation operator constructed in Section [6l Since we assumed A(s) > 
1/2 for all s, Lemma [6] implies that 

[H' s ,P]=0. (37) 

We can represent H' s as 

H' S = H + V, (38) 

where 

V = UlH Q U s -H + sUtVU s . (39) 

Lemma [2] implies that T> s has strength O(J) = 0(1). Applying Lemma U we conclude that 
sUsVU s has strength 0{J). Let us now focus on the term uIHqU s — Hq. We use 

U\H U S -H = -i f ds' U\, [V S ,,H ]U S ,. (40) 
J o 

Since T> s i has strength O(J), the commutator [T> s /,Hq] also has strength O(J). Applying 
Lemma Q] to the unitary evolution W,, we infer that uIHqU s — Hq has strength O(J). To 
conclude, we have shown that V' has strength O(J), and Eq. (I37p implies 

[y',p] = o, (4i) 

that is, V is a globally block-diagonal perturbation with strength 0{J). In Lemma [7] below, 
we will show that we can rewrite 

H' s = H + V = H + J2 X u (42) 
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where X u obeys [X u , P] = 0, where X u has strength O(J) and its support is near u. Applying 
Lemma [3] to Eq. (|42p implies that H' s can be written in the form 

H' a = H + V' = H + V" + A, (43) 

where V" is locally block-diagonal with strength O(J) and ||A|| decays faster than any power of 
L*. The statement of the theorem then follows straightforwardly from Lemma [5] which asserts 
that V" is relatively bounded by Hq with a small error decaying faster than any power of L* . 

The Lemma [7] that we need follows the idea in [17J to write a Hamiltonian of a gapped 
system as a sum of terms such that the ground states are eigenvectors of each term separately 
(in |16| a related idea of writing it so that the ground state was an approximate eigenvector of 
each term separately was considered). The properties of H' s that we use are that it is globally 
block diagonal, it has a spectral gap > 1/2, the perturbation V has strength J, and that it is 
unitarily related by U s to a Hamiltonian with a decay rate \i > 0. 

Lemma 7. Let H' s be defined as above. Then, we can re-write 

H' s = H + J2 X u (44) 

where X u obeys [X U ,P] = 0, where X u has strength O(J) and its support is near u. 

Proof. We start from representing V as V = J2ueA where V u includes only interactions 
affecting a site u. Then V u has strength J and its support is near u. We set 

/oo 
dtg(t) exp(iH' s t)V u exp(-iH' s t), (45) 
-oo 

where g{t) is a function satisfying g{—t) = g(t)* such that its Fourier transform g{oS) is infinitely 
differentiable, has <?(0) = 1, and g(uj) = for \oj\ > 1/2. Define 

/oo 
dt g(t) exp(iH' s t)Q A exp(-iH' s t). (46) 
-oo 

Then, 

/oo 
dtg(t)exp(iH>t)H' s exp(-iH' s t)= £ Qa + ^V u . (47) 

"°° Ag5(2) «GA 

By construction of g(co) we have (1 — P)QaP = (1 — P)V U P = 0. Hence both Qa and V u 
preserve P. Using the definition of V u we get 

/oo 
dtg{t) exp{iH s t)U s V u Ut exp(-iH s t). (48) 
-oo 

Recall that T> s has strength O(l). Thus we can apply Lemma Q] to the unitary evolution U s to 
infer that U s V u ul has strength O(J). Because g(oS) is infinitely differentiable, g(t) decays faster 
than any power. Also, by assumptions of the theorem, H s has strength 0(1) and decay rate 
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pi > 0. Hence we can apply Lemma [2] to the unitary evolution exp(iH s t) to infer that U s V u ut 
has strength O(J), Finally, applying Lemma Q] to the unitary evolution Jj\ we infer that V u 
has strength O(J). In addition, V u has support near u since all Hamiltonians obtained at the 
intermediate steps have support near u, see Lemmas 11 121 
We now consider the terms Qa- We have 

/oo 
dtg(t)exp(iH' s t)Q A exp(-iH' s t) (49) 
-CO 

/oo rt 
dtg{t) I dt 1 exp{iH' s t 1 )[V,Q A ]exp(-iH' s t 1 ) 
-oo JO 

/oo 
dbi f(h) exp (iH' s h)[V, Q A ] exp (-iflfo) 
-oo 

for some function f(t\) decaying faster than any power. It follows that 

/oo 
dh f(h) exp (iH s h)U s [V, Q A ]Ut exp {-iH s t x ). (50) 
-oo 

Applying the same arguments as in the analysis of Eq. ([15]) we conclude that Qa — Qa has 
strength O(J) and has support near u(A) — the center of the square A. In addition, Qa — Qa 
commutes with P since both Qa and Qa do. Let us define X u = V U + {Qa — Qa) where A is 
the square centered at u. □ 
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